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CN . Abstract. We examine the possibility to utilize in-medium charmonium for- 

mation in heavy ion interactions at collider energy as a probe of the properties 
of the medium. This is possible because the formation process involves recom- 
bination of charm quarks which imprints a signal on the resulting normalized 
I/"") ■ transverse momentum distribution containing information about the momen- 

tum distribution of the quarks. We have contrasted the transverse momentum 
spectra of J/ip, characterized by (pr 2 ), which result from the formation pro- 
ij cess in which the charm quark distributions are taken at opposite limits with 

regard to thermalization in the medium. The first uses charm quark distribu- 
qj ' tions unchanged from their initial production in a pQCD process, appropriate 

if their interaction with the medium is negligible. The second uses charm quark 
distributions which are in complete thermal equilibrium with the transversely 
expanding medium, appropriate if a very strong interaction between charm 
' quarks and medium exists. We find that the resulting (pr 2 ) of the formed J/ip 

should allow one to differentiate between these extremes, and that this differ- 
entiation is not sensitive to variations in the detailed dynamics of in-medium 
formation. We include a comparison of predictions of this model with prelimi- 
nary PHENIX measurements, which indicates compatibility with a substantial 
fraction of in-medium formation. 

Keywords: Quark-Gluon Plasma 
PACS: 25.75i.Nq, 12.38.-t 



1. Introduction 

The in-medium formation picture we consider here [[J uses competing formation 
and dissociation reactions in a Boltzmann equation to calculate the final J/ip pop- 
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ulation. The absolute value of this formation was found to be very sensitive to the 
underlying charm quark momentum distributions [ 2 . There is also quite strong 
variation of the results which depend on largely unconstrained model parameters 
involving details of the size and expansion profile of the deconfinement region. The 
initial PHENIX data suffered from low statistics, and was compatible with a 
fairly large region of model parameter space [0]. 

Recent work in this area concentrated on finding a signature for in-medium 
J I ip formation which is independent of the detailed dynamics and magnitude of the 
formation. We found that the px spectrum of the formed J/ip may provide such a 
signature 

The calculations which used initial charm quark momentum distributions from 
NLO pQCD amplitudes to generate a sample of cc pairs, were then supplemented 
by an initial-state transverse momentum kick to simulate confinement and nuclear 
effects. In the evolution of the interacting system size from pp to pA to AA colli- 
sions, the pt will be in general increased due to initial-state effects of interaction 
of constituents in the nuclei. 

{pt 2 )ab - (PT 2 )p P = A 2 [n A + n B - 2], (1) 

where ua is the impact-averaged number of inelastic interactions of the projectile 
nucleons in nucleus A, and A 2 is the square of the transverse momentum transfer 
per collision. The PHENIX measurements of J /ij) pt spectra in pp and minimum- 
bias d-Au interactions [ 6 allow us to determine the amount of initial state kx 
needed to supplement our collinear pQCD events. (This is equivalent to using a 
hadronization model similar to a color evaporation scheme, except that we assume 
that the resulting px of the resulting J/ip is determined by the pair px for all 
invariant mass of the combinations. One can then extrapolate to Au-Au and predict 
the spectrum of J /if) which are produced from hadronization of the initial " diagonal" 
cc pairs, again for minimum bias interactions. (We use diagonal to distinguish these 
pairs from the "off-diagonal" combinations which contribute to in-medium J /if) 
formation.) One finds ua = 5.4 from minimum bias d-Au interactions at RHIC 
energy (using <j pp = 42 mb), which leads to A 2 = 0.35 ± 0.14 GeV 2 . We note that 
the relatively large uncertainty comes entirely from the difference in px broadening 
in d-Au between positive and negative rapidity. 

Our prediction for the "normal" evolution of the px spectrum in Au-Au inter- 
actions is shown by the triangular points in Fig^ 

The properties of the normalized pt spectrum for the formation process follow 
from two separate effects: First, the fact that the process is dominated by the off- 
diagonal pairs introduces a modified initial px distribution. Next, one weights these 
pairs by a formation probability for J/ip- We use the operator-product motivated 
cross section for cc forming J /tp with emission of a final-state gluon, which of course 
is just the inverse of the dissociation process. However, any cross section which has 
the same general properties as this one gives essentially the same result [ |o] . We 
show by the square points in Fig^the prediction for the formed J/ip- One sees that 
this spectrum is substantially narrower than the one with no in-medium formation. 
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Fig. 1. Comparison of in-medium J ftp transverse momentum spectra predictions 
from various scenarios. 



The same procedure was employed when using quark momentum distributions 
which follow for charm in thermal equilibrium with the expanding region of de- 
confinement. The parameters of temperature and maximum transverse expansion 
rapidity are determined by a fit to this thermal behavior of the produced light 
hadrons. The application to charm quarks was originally motivated in Ref. [Cj, 
who showed that the low- px spectrum of decay leptons from charmed hadron de- 
cays would not be able to differentiate between the thermal and a purely pQCD 
distribution. We show here, however, that the px spectrum of in-medium formed 
J/ip is very sensitive to this distribution. The circles in Fig^ result from formation 
calculations using T = 170 MeV and yTmax = 0.5 for the thermal charm quarks. 
One sees that this pt spectrum is narrower yet than in-medium formation from 
pQCD quarks. Finally, we show by the stars the pr spectrum of J/iJj which them- 
selves obey this thermal distribution. The resulting spectrum falls between the 
in-medium formation spectra for either pQCD or thermal charm quark momentum 
distributions. 



2. Centrality behavior 



We now proceed to investigate the variation of the pQCD-based results with respect 
to the collision centrality in Au-Au interactions. First, we use the value of A 2 
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extracted from pp and pA data, together with values of t%a calculated as a function 
of collision centrality, to recalculate the (pr 2 ) values for either the initial production 
or the in-medium formation separately. This provides the centrality behavior of 
the J/ip spectrum in the case that one or the other of these mechanisms is solely 
responsible for the total J/ip population. We show these results together in Fig. 
13 One sees as expected that (pr 2 ) is maximum for the most central collisions, but 
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Fig. 2. Centrality dependence of J /4>{pt 2 ) contrasting predictions assuming either 
100 % production from initial cc pairs or 100 % in-medium formation. 

the absolute magnitudes are widely separated for initial production and in-medium 
formation at each centrality. One should note that the uncertainties are dominated 
by the difference between the px -broadening measurements at positive and negative 
rapidities in the d-Au interactions. Thus the point-to-point uncertainties are much 
smaller for the centrality behavior. We have also included separate values for (pr 2 ) 
in the region limited by a maximum px of 5 GeV, to facilitate comparison with 
experiment in this same range. 

There exist preliminary results from PHENIX for (px 2 ) of J/tp produced in Au- 
Au collisions [|H] ■ These are reported as a function of the number of nucleon-nucleon 
collisions N co u. In order to compare with our predictions, we transform centrality 
to Ncou using a Glauber model. The resulting predictions are shown in Fig. EH 
Although there is substantial uncertainty in the absolute values, one can infer a 
clear preference for in-medium formation over the initial production prediction. 
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Fig. 3. Comparison of (pr 2 ) predictions for initially-produced J/ip and in-medium 
formation with preliminary PHENIX measurements. The collision centrality is pa- 
rameterized by the number of nucleon-nucleon collisions N co u. The theory curves 
and uncertainties are taken from Fig. [21 

In order to provide a meaningful prediction for the overall J/ip spectrum, one 
should of course include both initial production and in-medium formation together 
as sources. This requires some estimate of the relative magnitudes of these processes, 
and is subject to considerable model uncertainties. What we can say, however, is 
that in-medium formation will be most dominant for central collisions, where the 
quadratic dependence on N c5 is enhanced. Conversely, one expects that initial 
production will increase in relative importance for very peripheral collisions. To 
get an approximate idea of how this effect will appear, we revert to our original 
model calculations which included the absolute magnitude results One relevant 
parameter is the number of initial pairs, N c£ , parameterized by its value at zero 
impact parameter. These results are shown in Fig. 0]for three representative values 
of N c5 — 10, 20, and 40, which span the range of values extracted from STAR 
and PHENIX [ measurements of charm inferred from semileptonic decays and 
direct reconstruction of D mesons. 

The 100 % curves for initial production and in-medium formation remain as in 
Fig. |21 Between these two extremes are shown the calculations including centrality- 
dependent fractions of each, as described above. Each N c5 value is shown by the 
solid, dashed, and dot-dashed lines, respectively. The duplication of these lines 
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Fig. 4. Variation of predictions in Fig. [3] resulting from a model calculation of 
the relative fraction of J/ip which originate from initial production or in-medium 
formation. 

results from different initial temperature values of 300, 400, and 500 MeV, which 
provides variation in the gluonic dissociation rates and the deconfinement lifetime 
in the model calculations. One sees that, with the exception of extremely peripheral 
collisions, all of the predictions remain substantially below that for initial formation 
alone. It would require a substantial reduction of experimental uncertainty in the 
measured (pt 2 ), however, to pin down the preferred model parameters appropriate 
for in-medium formation. 

3. Hadronization Contribution 

In addition to in- medium formation of J/ip, one needs to consider the possibility 
of subsequent production at the hadronization transition. This contribution must 
be significant, as can be confirmed by an examination of the absolute values of 
in-medium formation in the model calculations. For central collisions, our calcu- 
lated values of the ratio Nj/^/N cS typically range from 0.004 to 0.01 as the model 
parameters are varied. This means that essentially all of the initial cc pairs will 
survive through the deconfined phase and some fraction of these must hadronize 
into additional J/ip. For an estimate of this fraction, we utilize the statistical 
hadronization model [ • The most recent applications of this model to RHIC 
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Au-Au collisions [ 1121 IT3| predict values for Nj/,p/N cS m 0.005, comparable to that 
for our in-medium formation. Thus we need to consider contributions to the final 
J I ip population which originate in approximately equal amounts from each mecha- 
nism. Fortunately one can make this combination for central collisions alone, since 
both in-medium formation and statistical hadronization predict the same centrality 
behavior [EJ. The results of such a calculation are shown in Fig. 03 
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Fig. 5. Comparison of px spectra for J jtp resulting from variable fractions of 
in-medium formation and subsequent hadronization of residual charm. 



Shown are the px spectra for the in-medium formation utilizing pQCD charm 
distributions (diamonds) and thermal plus flow charm distributions (squares). They 
are at the extremes of (px 2 ), as previously noted. Also shown is the thermal plus 
flow distribution for a direct J ftp, which we take as our expectation if statistical 
hadronization is the dominant source of J/tp- The distributions which fall between 
the extremes are the result of various (but comparable) contributions of in- medium 
formation and hadronization. One sees that these combinations predict the spectra 
evolve toward that for the thermal J/tp as its fraction increases. However, the 
rate of this evolution is substantially greater for in-medium formation from thermal 
charm quarks than the corresponding behavior for in-medium formation from pQCD 
charm quarks. Thus the possibility to utilize the px spectra as a probe of charm 
thermalization in the medium remains a viable option. 
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